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ABSTRACT 

In this work we present a semi-analytical model of galactic wind for dwarf 
galaxies where thermal and turbulent /momemtum driving mechanisms are 
studied. The model takes into account wave and internal adiabatic heating 
mechanisms, as well as radiative and adiabatic cooling. The importance of 
external sources of energy is discussed. We also studied the role of the spatial 
distribution of dark matter in the acceleration of the wind and on the mass 
loss rates. The basic model parameters that regulate the wind mass loss rate 
and terminal velocity are the gravitational potential profile, the equilibrium 
temperature of the gas and the amplitude of the turbulent motions of the gas. 
We found that dwarf galaxies are likely to present quasi-stationary winds with 
mass loss rates larger than lCP 3 M Q yr -1 even in the absence of turbulent 
motions (which is possibly related to the SNe feedback) , if the interstellar gas 
is heated to T > 10 4 — 10 5 K. We also found that the wind mass loss rate is 
enhanced for cusped dark matter distributions, such as the NFW-profile, due 
to the increased pressure gradients at the center of the galaxy. The solutions 
presented here may serve as benchmarks for numerical simulations, and as 
inputs for single zone chemical evolution models of dwarf galaxies. 
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1 INTRODUCTION 

Outflows at galactic scales are ubiquitous in the Universe, typically related to the most 
energetic physical phenomena known, such as accretion into active galaxy nuclei (AGNs) 
and supernovae (SNe) driven winds (Veilleux, Cecil & Bland-Hawthorn 2005, Chen et al. 
2010, Steidel et al. 2010, Erb et al. 2012, Martin et al. 2012). Much less magnificent (and 
energetic) than these, but not less relevant, winds from less massive galaxies also play an 
important role on the chemical evolution of the intergalactic medium (IGM), and on the 
physical properties of the interstellar medium (ISM) of the galaxy itself (Dekel & Silk 1986; 
Bradamante, Matteucci & D'Ercole 1998; Fragile et al. 2003; Tremonti et al. 2004; Kirby, 
Martin & Finlator 2011; Peeples & Shankar 2011, and others). 

Few tens of dwarf elliptical and spheroidal galaxies have been observed in the Local 
Group, most of them orbiting the Milky Way and M31. These are small objects, typically 
< lkpc, but relatively massive, presenting dynamical masses of M* ~ 10 7 — 1O 9 M , being a 
small fraction visible M* ~ 10 5 — 10 7 M Q (Mateo 1998). The large mass-to-luminosity rela- 
tion obtained indicates that these objects are dominated by dark matter (e.g. Armandroff, 
Olszewski & Pryor 1995). 

Cosmological numerical simulations predict the existence of hundreds to thousands of 
low mass dark matter haloes orbiting large galaxies at groups as ours. In a possible scenario, 
most dark matter haloes may be associated to extremely low luminosity stellar components, 
which are simply not detected by current instruments (Somerville & Primack 1999). In con- 
trast to this, some numerical simulations predict a large number of merging and disruption 
of sattelites in the Local Group (e.g. Bullock & Johnston 2005, Robertson et al. 2005). 
Galactic winds may remove the gas from dwarf galaxy haloes quenching star formation. The 
systematic low metallicity observed in dwarf galaxies also indicates that the star formation 
process is quenched at a given point of the galaxy evolution (Kirby, Martin & Finlator 2011), 
possibly related to winds. 

In the vast majority of galactic wind models star formation and SNe are assumed as main 
sources of energy and momentum (see e.g. Larson 1974; Chevalier & Clegg 1985; Bradamante, 
Matteucci & D'Ercole 1998; Stringer et al. 2012, Ruiz et al. 2013). Basically, these models 
assume that the SNe energy is thermalized at the ISM, increasing the gas thermal pressure 
that drives the galactic wind. The energy released by the SNe, or a fraction of it (typically 
~ 1 — 10%; Silich, Tenorio-Tagle & Munoz-Tunon [2003]), may then be directly compared to 



Galactic winds from Dwarf Galaxies 3 

the binding energy of the gravitational potential (e.g. Bradamante, Matteucci & D'Ercole 
1998; Lanfranchi & Matteucci 2007, Cote et al. 2012). This method is very simplified and 
do not account for different spatial distribution of the gas and the star forming regions, 
temperature distributions, and processes that may influence the dynamics of the wind, such 
as local turbulence and radiative losses. 

A complete description of the dynamics of SN ejecta interacting with the local ISM, 
and further acceleration of the galactic wind, is a very complex problem. Heating from 
cosmic rays, fractional ionization and radiative cooling are only few of the processes that 
make the effort of modeling the galactic winds even more difficult. Numerical simulations 
are able to describe many of these processes occurring in galactic outflows simultaneously 
(e.g. Recchi, Matteucci & D'Ercole 2001; Mori, Ferrara, & Madau 2002; Fragile et al. 2003; 
Fujita et al. 2004; Ruiz et al. 2012). However, different schemes and numerical resolutions 
provide different results. It is well known that SPH schemes are unable to evolve shocks 
and certain instabilities, such as the Rayleigh- Taylor and Kelvin-Helmotz, properly. Even 
numerical recipes for injecting SNe energy in simulated galaxies are still under debate, from 
the randomness of starbursts in time and location, to the actual efficiency on the conversion 
of the SN thermal to ISM kinetic energy. For example, Recchi et al. (2001) simulate a single 
starburst with all stars at a narrow location, as well as Fujita et al. (2004), and assume that 
3% of the type II SNe energy feedback would thermalize with the ISM, which results in a 
single and gigantic hot bubble that drags all gas out of the galaxy. Mori et al. (2002), on the 
other hand, assumed a net net injection of energy into the ISM of ~ 30% of the SNe energy, 
but obtained similar results. Ruiz et al. (2012), on the other hand, showed that the buoyant 
instabilities may be responsible to up to 40% of the mass loss of a dwarf galaxy. It is crucial 
therefore to provide the numerical community analytical (or semi-analytical) benchmarks. 

In this work we present an analytical model of galactic winds in order to study the 
role of different dark matter density distributions, temperature distribution, and turbulent 
pressures, on the gas acceleration. The model is applied for dwarf spheroidal galaxies. We 
particularly focus on the possibility of having strong winds, i.e. with mass loss timescales 
shorter than lGyr, in the absence of starbursts, for which thermally driven winds would 
suffice. We also study the scaling of the mass loss rate of the hybrid thermally /wave driven 
winds with the energy density of the kinetic perturbations, which may be correlated to star 
formation activity. The model is described in Section 2. In Section 3 we present the results 
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where we compare different scenarios, followed by a brief discussion of their impact on the 
evolution of dwarf galaxies, in Section 4, and by the Conclusions. 



2 BASIC EQUATIONS 

If rotation does not play relevant roles in the dynamics of the interstellar gas we may assume 
radial symmetry for this problem, i.e. 8$ — — 0. The basic hydrodynamic equations to 
describe the gas are the mass continuity and momentum equations, given by: 

^ + V-( P u) = 0, (1) 
du „ GM(r) ldP 

— + u-Vu = — + g rad + g wav , (2) 

at r 2 p or 

where u, P and p represent the plasma velocity, thermal pressure and mass density, respec- 
tively. The additional terms g ra d and g wav correspond to the specific net force generated by 
the radiation field and by waves, respectively. In this work we will also neglect the role 
of radiation pressure in the acceleration of the galactic wind, but will take into 
account the radiative losses in the energy equation. 

Typically, models of galactic winds assume the feedback from SNe as the main 
driver of such kinetic power. The injection of kinetic turbulence instead of the 
dynamical evolution of energy blasts when simulating SNe feedback have been 
discussed, e.g. by Falceta-Gongalves et al. (2010) and Scannapieco & Bruggen 
(2010). Basically, the advantage of their recipe for SNe feedback is to bypass the 
thermal evolution of the SNe Sedov's phase. After this quasi-adiabatic expansion 
phase most of the energy is lost due to radiation, but part of the initial energy 
has been also converted into kinetic. The stochastic distribution of several SNe 
and star formation regions therefore would result in a turbulent medium. How- 
ever, this prescription requires the adoption of an efficiency factor (rj) for the 
conversion of the supernova energy into turbulence that is not easily estimated. 
As pointed by Ruiz et al. (2013), the generation of SN driven turbulence de- 
pends on the distribution of the blast waves in time and space. This because 
the "snow plow" phase stalls at scales of order of 10 - 100 pc, depending on the 
local ISM gas density and temperature. As an example, shown by their numeri- 
cal simulations, the merging of hot cavities efficiently generate Rayleigh-Taylon 
unstable regions, while SNe localized far from the central region may brake the 
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ISM and only a minor fraction of its energy is thermalized with the ISM. These 
effects result in a different scenario as the result of simultaneously and centrally 
localized packs of SNe. 

In this paper, when solving the momentum and energy hydro dynamical equa- 
tions, we will deal with the wave term as a kinetic pressure only. We do not at- 
tempt to relate these waves a priori to any specific physical mechanism occuring 
in the galaxy, such as turbulence generated by gas infall or tidal forces, or by 
any internal source, such as supernovae and stellar winds. 

The acceleration provided by the wave pressure gradient is given by: 

where e represent the energy density of the waves, and is related to the averaged squared 
amplitude of the acoustic perturbations ((SV 2 )) as: 

e = p(5V 2 ). (4) 



The energy equation, given by: 

d (u 2 5k B T TT \ ude e(r) . , . , _, 

»"ai [T + — ~ u " ) + 2Tr = l + c *> - + r - < 5 > 

where A(T) is the cooling function, T the net heating source from external sources, c s the 
sound speed, U grav the gravitational potential and L(r) the damping length of the acoustic 
waves, as a function of the radius. 

The set of equations is closed by the differential equation that describes the gradient of 
the wave energy density. Under the WKB approximation, the energy density of the waves 
follows the relation: 

where we have defined the sonic Mach number M. = u/c s , and the index "0" represents the 
variable value at r = r . Therefore, we may now obtain the differential form: 
de e e(l + 3M) fldu 1 dT\ 

Since we are mostly interested in the stationary solution of the hydrodynamic equations 
we will also consider d t = 0. Now substituting Eqs. 1 and 7 in the momentum and energy 
equations (Eqs. 2 and 5), we obtain: 
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u dr 



(9) 



respectively, with Q w = e(u + c s )/L the total heating from the decaying energy of the waves. 

The set of equations above can now be integrated spatially in order to give the radial 
distributions of the wind physical parameters. 



3 RESULTS 

The physical processes that bring much complexity to the modeling of the galactic winds are 
related to the radiative heating and cooling of the plasma. Cosmic rays, ionizing radiation 
and the photoelectric effect on dust particles represent the main heat sources of the Galactic 
interstellar gas. These processes however are stricktly dependent on the stellar population 
itself, as well as on the dust properties (for the photoelectric effect) and magnetic field (for 
the diffusion of cosmic rays) of the ISM, making any attempt to quantify their contributions 
in early dwarf galaxies not practical analytically. 

In this sense, we will consider in this work two special cases. Firstly, we will assume the 
plasma to be isothermal. This is equivalent to assume that the timescales for radiative cooling 
and heating are much shorter than those of the wave heating and adiabatic cooling due to the 
expansion, and that the equilibrium between both radiative cooling and heating occurs at 
a given equilibrium temperature T cq . Secondly, we consider the full non-isothermal model, 
but assume that the external heating sources are irrelevant compared to the other terms 
of Equation 9, i.e. r — > 0. As discussed later, this is particularly interesting because even 
though not modeled properly, the heating sources (r) are shown to be extremely important 
for sustaining a steady galactic wind. 
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3.1 The special case of isothermal galactic winds 

In order of analytical simplicity let us start with the simplest case with no temperature gra- 
dients (dT/dr = 0), i.e. an isothermal galactic envelope. This assumption is not unrealistic 
in two conditions, when the strong UV radiation field of a starburst (or cosmic rays) is able 
to ionize most of the volume of the ISM of a dwarf galaxy, except for the clustered molecular 
clouds where star formation will continue to occur, or if the galaxy has no stars formed yet 
and the gas is basically molecular. In this case, the momentum equation (Eq. 8) is written 
in the form: 



ldu 

u dr 



u 



2 (5V 2 ) 1 + 3M 
° s 2 1 + M 
GM(r) 2c| {5 V 



2L 



(10) 



where the isothermal sound speed c s = (ksT/ fimu) 1 ^ 2 is used. 



3.1.1 fore^O 

Let us first consider the case where the energy density of acoustic waves is negligible com- 
pared to the thermal and gravitational energy densities. In this case we reduce the galactic 
wind problem to the standard Parker solution for a stellar wind: 

(2c 2 s GM(r) \ 

ldu = {~ —) ^ 

u dr (u 2 — c 2 s ) 

Typically, for stars, the numerator at the fraction in the right of Eq. 11 is negative at the 
lower boundary of the wind. Thus this equation has an unique solution when both numerator 
and denominator of the right side fraction must vanish, at the critical radius. This position 
is the so-called sonic point, since w(r crit ) = c s . The critical radius is obtained by making the 
numerator vanish. 

Now, considering a galaxy instead, the position of r crit depends on the mass of the 
galaxy M(r cr ; t ). The dark matter component dominates the gravitational potential of dwarf 
galaxies. 

The mass distribution of the dark matter halo may be described by a Navarro-Frenk- 
White (NFW) profile: 

M NFW (r) = AirpoRl 



.V Rft J i^ + rJ' 
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Figure 1. Minimal core mass density of the dark matter component needed for an accelerating wind at r = r cr jt l, as a 
function of the characteristic radius Rg, obtained for T = 10 4 K. The relation has been calculated for r cr ; t j = lpc (solid), lOpc 
(dash-dot), 50pc (dashed) and lOOpc (dotted). 

where Rs represents a characteristic radius and po the core mass density of the dark matter. 
M NFW diverges as r — > oo so we have to truncate the NFW profile at a given radius -Rtmno 
which is assumed as lkpc for this work. For larger radii the total enclosed mass of the galaxy 
will be constant, i.e. M NFW (i? trunc ). For the sake of simplicity we will only consider the dark 
matter mass contribution to the gravitational potential. This approximation is valid for all 
the cases we study in this paper, as the ratio of the dark matter to the gas potential is 
much larger than unity, for the radii computed (r > lOpc). This approximation is important 
also because once the galactic wind is set the gas mass profile becomes time-dependent, 
in contradiction to the steady-wind approximation used in this work. Similarly, we are not 
considering cosmological evolution of the halo dark matter mass or characteristic radius. 
Therefore we do not vary the characteristic radius, profile or total mass of the dark matter, 
function of redshift. 

At small radii m<c s and the gas is accelerated outwards only if the numerator of this 
equation is negative. A major issue to the theory of pure isothermal winds in dwarf galaxies 
arises from the fact that the gravitational potential may be small at the central regions of 
the galaxy. As a consequence, the numerator of Eq. 11 will be positive resulting in d r u < 0, 
i.e. the wind is actually a decelerating breeze. This will occur up to larger distances where 
the acceleration of gravity becomes dominant compared to the thermal component, and d r u 
becomes positive. We define the position where this transition occurs as the first critical 
point (r critj i). 




— r crit.1 = 1P C 



— 10pc 

- ■ 50pc 
lOOpc 
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For a NFW dark matter distribution the critical location of the wind is obtained by 
solving the following equation, obtained by setting the numerator of Eq. 11 equals to zero: 



where A = 27iGp R s c~ 2 and x = r CTit /R s . This equation may have up to 2 real solutions, 
as we discuss below. Depending on the gas temperature, if the core density po is lower than 
a given threshold there is or only 1 root for the above equation, which means that the gas 
is not thermally accelerated. The galactic wind will be stationary only in the cases where 
there are 2 real roots. 

The first critical point depends on T, p and R S - In Figure 1 we show the dependence of 
the minimum density p required, as a function of Rs, for the wind to be accelerated from 
^crit.i, for T = 10 4 K. For example, if the dark matter distribution has p = lCT 21 g cm~ 3 and 
Rs = 300pc, the galactic wind starts to be accelerated at r ~ 2pc. Basically, at r > r crit l 
the thermal energy is converted into kinetic as long as the temperature is kept constant. 

Where the wind velocity reaches the sound speed the denominator of Eq. 11 vanishes, 
and a real solution for this problem is only obtained if the numerator vanishes as well. We 
define this position as the second critical point for the isothermal galactic wind r critj 2- This 
location represents the second root of Eq. 13. 

The full solution of u(r) is obtained by integrating Eq. 11 numerically, for which we 
used a 4 th order Runge-Kutta integration scheme. The calculations are performed from 
r = lOpc, with a spatial resolution of 10~ 2 pc, up to r max = 3000pc. The mathematical 
solutions are independent on the minimum and maximum radii set, given that the sonic 
points are present in the chosen range. However, since energy injection should occur mostly 
at scales smaller than few tens of parsecs (e.g. starbursts) it is actually irrelevant to study 
the thermal acceleration at the most central regions of the galaxy. The gas temperature, 
Rs and p are initially set as constants. In order to obtain the critical solution at r criti2 we 
varied the initial condition uo- 

Two different models are presented in Figure 2 to illustrate the acceleration behavior of 
the thermal galactic wind. For these we set Rs = 200pc for both models, but T = 10 4 K and 
Po = l(T 22 g cm' 3 (solid) and T = 10 5 K and p = l(T 21 g cm" 3 (dashed). We plot both the 
sonic Mach number (u/c s ) and the gradient of wind velocity d r u as a function of the radius. 

The galaxy is actually surrounded by the intergalactic medium, which eventually becomes 
dynamically important for the expansion of the wind. The stagnation of the wind occurs 




(13) 
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Figure 2. Wind sonic Mach number M s = u/c s and gradient of wind velocity d r u as a function of the radius, calculated for 
two models with T = 10 4 K and p = 10 _22 g cm -3 (solid), and T = 10 5 K and p = 10 _21 g cm -3 (dashed). R s = 200pc for 
both cases. The third model shown represents the wind Mach number for a logarithmic profile of the dark matter density (see 
text). 

where the kinetic pressure of the wind is approximately balanced by the thermal pressure 
of the intergalactic medium, i.e. p wm dw;L> ~ ^igm&^igm- F° r ^igm = 10~ 3 cm~ 3 , T = 10 6 K 
and Mqo = 10 — 30km s— 1, we obtain p w i n d ~ 10~ 25 — 10~ 26 g cm -3 , which corresponds to 
^stag ~ 500 — lOOOpc. For the local Universe, such large IGM gas densities may occur at 
cluster of galaxies and extended gas haloes of massive galaxies, which the dwarfs may orbit. 
At high redshifts (e.g. z > 5), deep gravitational potentials, such as clusters, do not exist 
yet and massive galaxies are probably just forming. In this case the wind of the low mass 
objects will basically interact with the diffuse IGM. The average density of baryons evolves 
with the scale factor a as n(a) SH^^STrGrripa 3 )^ 1 ~ 1.1 x 10 _5 (1 + z) 3 Q b h 2 cm~ 3 , and 
^igm = 10~ 3 cm~ 3 may be obtained for z ~ 9. The asymptotic wind velocity obtained is 
typically few times larger than the sound speed, depending on the model. Therefore if the 
galactic wind stagnates at 500 < r < lOOOpc because of the pressure of the intergalactic 
medium (IGM) the wind velocity is not able to reach maximum values, and Uoo ~ c s becomes 
a good approximation. 

The gas mass loss rate is obtained by, 
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M = Anr 2 Pgas (r)u(r), (14) 

which may be calculated at any position of the wind due to mass conservation, e.g. at 
r = 200pc, where the models shown in Figure 2 result in M ~ 3.9 x 10~ 5 n2oo(cm -3 ) and 
1.7 x 10~ 3 x -n,2oo(cm~ 3 ) M Q yr _1 for T = 10 4 and 10 5 K, respectively. As a result, after 
~ lGyr, the total gas mass lost by an isothermal galactic wind of a typical dSplJl] would be 
of ~ 0.5% and 20% for the models with T = 10 4 and 10 5 K, respectively. 

3.1.2 The role of core or cusped dark matter profiles 

Recent observational and numerical works suggest that the NFW profile may not be the best 
description for the distribution of dark matter in dwarf galaxies (e.g. Burkert 1995; van den 
Bosch et al. 2000; de Blok & Bosma 2002; Kleyna et al. 2003; Simon et al. 2005; Walker et 
al. 2009; Governato et al. 2010; Oh et al. 2011; Del Popolo 2012; Jardel & Gebhardt 2012). 
Logarithmic profiles have been shown to match better the stellar velocity dispersion of the 
dSphs of the Local Group (see Persic, Salucci & Stel 1996, Salucci et al. 2012). A different 
dark matter distribution can result in different accelerations of the galactic wind. In order 
to test this possible effect, and to compare to the standard NFW profile, we also consider a 
logarithmic gravitational potential for the dark matter density distribution which results in 
the following density and total enclosed mass distributions: 

'-M-i^r^!) 5 ' < 15 > 

V 2 r 3 

M los (r) ~ -, (16) 

G r^ + r z 

where V s represents the virial circular orbital velocity at r s , the characteristic radius for the 
dark matter distribution. 

The advantage of using Eq. 16 is that the critical radii are obtained analytically. The 
first critical radius is given by: 



crit I V} - 2c? 



o r 2 r 2 \ a 

" ~ • 1 ~ 45 pc, (17) 



considering V s = 100km s -1 , r s = lOOpc and T = 10 5 K. 

It is easy to notice that there is no second critical radius for the logarithmic profile. 
However, a saturation of the dark matter mass is expected to occur somewhere at r > r s , 
mostly due to tidal disruption within the group of galaxies. If we define a cutoff radius for 

assuming an initial gas mass of 10 7 M Q . 
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the dark matter distribution at r = 300 — 500pc then the second critical radius occurs at 
r crit ,i ~ 1.4 - 2.9kpc. 

For the logarithmic dark matter density profile the wind terminal velocity is similar to 
the one obtained with a NFW profile, as shown in Figure 2. For T = 10 5 K, we obtained 
asymptotic terminal velocities w^f w ~ 3.8c s while u^ g ~ 3.2c s , for a cutoff at r = 500pc. 
The mass loss rate, on the other hand, is considerably different. We obtained M ~ 2.1 x 
10 _6 n2oo(cm -3 ) M Q yr _1 , which is 3 orders of magnitude lower than the obtained with a 
cusped NFW profile. The reason for this difference is that the wind is accelerated at larger 
radii for a log-profile, where the gas density is much lower. 

3.1.3 The case of e ^ and L ->■ oo 

SNe, stellar winds, tidal stirring and internal dynamic motions of the galaxy may drive 
turbulence in the interstellar medium. Once these turbulent waves are present the radial 
gradient of their kinetic pressure accounts for an extra acceleration term. This gradient is 
strongly dependent on the wave damping. If we assume that the waves propagate with no 
specific process of energy loss, except for the spatial dilution of the waves as they propagate 
outwards, i.e. L — > oo, Equation 10 results in critical points which are exactly the same 
as in the case with no waves. However since the wind velocities now must be larger at the 
second critical point because of the extra term (SV 2 ), an increase in the mass loss rates and 
terminal wind velocities are expected. 

In Figure 3 we present the wind velocity profiles considering T = 10 5 K and a NFW 
profile with p = 10~ 21 g cm -3 (dashed) and R s = 200pc, for wave amplitudes (5V 2 } = 
(dotted), 0.1c 2 (dot-dashed), 0.2c 2 (dashed) and 0.5c 2 (solid). The wind terminal speeds and 
mass loss rates are also shown in Table 1. 

As expected, the wind mass loss rate increases with the initial energy density of the 
waves. As shown in Table 1, the mass loss rate is basically proportional to the turbulent 
Mach number for (5V 2 ) /c 2 > 0.1. A linear regression analysis provide M oc (5V 2 )^' 88 , 
with reduced x 2 = 0.96. The increase in the mass loss rate with the wave energy density is 
noticeable in both dark matter profiles. For (SV 2 ) /c 2 s < 0.1 the role of acoustic waves is 
irrelevant. 

The extra acceleration of the wind from waves below the first critical point results in 
increased mass loss rates. In order to continue presenting a physical stationary solution 
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Table 1. Description of the models for isothermal winds 



Model 


po(g cm 3 ) 


-Rs(pc) 


V s (km s" 1 ) 


r s (pc) 


T(K) 


4 


L(pc) 


c s 


"200 v u 


1 3\ 
'cm ) 


IsoNFWl 


10 


22 


200 


- 


- 


10 4 


0.0 


CO 


2.1 


3.9 x 10" 


-5 


IsoNFW2 


10~ 


21 


200 


- 


- 


10 6 


0.0 


CO 


3.S 


1.7 x 10" 


-3 


IsoNFW3 


10" 


21 


200 


- 


- 


10 5 


0.1 


CO 


4.2 


2.2 x 10" 


-3 


IsoNFW4 




21 


200 


- 


- 


10 6 


0.2 


CO 


4.1 


3.0 x 10" 


-3 


IsoNFW5 


10- 


21 


200 


- 


- 


10 5 


0.5 


CO 


4.3 


6.4 x 10" 


-3 


isoin r vv u 


10" 


21 


9nn 
zuu 






1 n"' 


1 n 
1 .u 




4 1 


1 7 y 1 fl" 


-2 


IsoNFW7 


10" 


21 


200 


- 


- 


10 5 


0.5 


20 


3.4 


1.6 x 10" 


-3 


IsoNFW8 


io- 


21 


200 


- 


- 


10 5 


0.5 


50 


3.5 


1.7 x 10" 


-3 


IsoNFW7 


io- 


21 


200 






10 5 


10.0 


20 


3.4 


1.6 x 10" 


-3 


IsoNFW8 


io- 


21 


200 






10 s 


10.0 


50 


3.6 


1.9 x 10" 


-3 


IsoNFW9 


io- 


21 


200 






10 s 


10.0 


100 


3.5 


6.1 x 10" 


-3 


IsoLogl 








100 


100 


10 5 


0.0 


CO 


3.2 


2.1 x 10" 


-6 


IsoLog2 








100 


200 


10 s 


0.0 


CO 


2.4 


6.0 x 10" 


-6 


IsoLog3 








100 


200 


10 5 


0.2 


CO 


2.6 


1.5 x 10" 


-5 


IsoLog4 








100 


200 


10 s 


0.5 


CO 


2.7 


5.9 x 10" 


-5 


IsoLog5 








100 


200 


10 s 


1.0 


CO 


2.8 


6.4 x 10" 


-4 


IsoLog6 








100 


200 


10 5 


1.0 


20 


3.2 


7.4 x 10" 


-4 


IsoLog7 








100 


200 


10 5 


10.0 


20 


3.3 


7.1 x 10" 


-4 


IsoLog8 








100 


200 


10 5 


10.0 


100 


3.3 


2.5 x 10" 


-3 



1E+01 
1E+00 
1E-01 
1E-02 
3 1E-03 
1E-04 



1E-05 



1E-06 
1E-07 

10 100 1000 

r (pc) 

Figure 3. Wind sonic Mach number A4 S = u/c s as a function of the radius, calculated for models with wave amplitudes 
(SV 2 ) Q = (dotted), O.lc^ (dot-dashed), 0.2c^ (dashed) and 0.5c^ (solid). We adopted T = 10 5 K, and a NFW profile with 
p = 10~ 21 g cm~ 3 (dashed) and R s = 200pc 
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at the second critical position the initial velocity of the obtained solution increases with 
(SV 2 ), which results in increased mass loss rates. In a marginally supersonic regime, i.e. 
(SV 2 y ,z ~ c„ the mass loss rate is increased by an order of magnitude. On the other hand, 
the waves play a minor role on the terminal wind velocities. This because the second critical 
point is located too far from the central region of the galaxy and, at that location, the energy 
density of the waves is too low. It means that beyond the second critical position, where 
further acceleration (above a value very close to the sound speed) should occur, the waves 
have little energy to effectively work over the wind. 



3.1.4 The case of e 7^ and finite L 

The picture described above changes when we allow the waves to decay. As the waves are 
damped the radial gradient of the wave energy density becomes steeper. As a consequence, 
the net force generated by the wave kinetic pressure over the plasma increases. However, 
a not so obvious result is that this extra pressure actually helps to radially decrease the 
wind speed if r < r crit l . Within this region, as discussed above, d r u is negative and a short 
damping length L worsens this scenario (see Eq.10) by reducing u even more. Typically, 
the wave damping length is short enough for the wave energy to be mostly released at radii 
smaller than that of the acceleration region, i.e. L < r crit l . At r > r crit 2 the wind stationary 
solution becomes similar to the standard isothermal case, as given in Table 1. 

Actually, our results show that the role of the wave damping depends on the dark matter 
profile that is considered. For a NFW profile the mass loss rates for the models without waves 
are very similar to those with damped waves, for L < 50pc. If L > 50pc the mass loss rates 
of the models with damped waves become similar to those without damping. This relation 
is not obtained for the log-profile though, for which the models with and without wave 
damping present similar mass loss rates and wind terminal velocities. 

Obviously, the model of an isothermal wind with no wave pressure is a crude approxima- 
tion of a much more complex phenomenon. Nevertheless, this model represents an interesting 
analytical benchmark for a more complete understanding of the problem of mass loss in dwarf 
galaxies. The full solution of the thermal wind as given above, i.e. non-isothermal and in- 
cluding the net forces of acoustic waves within the galaxy, are presented in the next section 
where we compare those results to the simplest model described above. 
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3.2 Non-isothermal winds 

Now let us assume that the temperature gradient is not zero. Unfortunately, it is not possible 
to consistently describe all the external heating sources that may occur at the early stages 
of dwarf galaxies. Therefore we will simply assume that T in Eq.9 is irrelevant compared 
to the other terms (at least for T > 10 4 K), namely the radiative cooling, wave heating and 
adiabatic expansion. This is possibly not a good approximation from the time when a large 
number of stars is already formed because their luminosity may become dominant in the 
energy budget of the interstellar gas. Eqs. 8 and 9 combined now lead to: 



1 du 
u dr 



2 5cj (5V 2 ) 1 + 3M 
"3 2 1 • M 



GM(r) + 2cf 



4 (SV 



2\ 



+ 7T" + 



3r L 



1 2/.M + 1 



2 3 V M 



The cooling function A(T) used in our calculations was obtained from the interpolation 
of a cooling efficiency table for an optically thin gas for T ^ 10 4 K (Gnat & Sternberg 2007), 
for a low metallicity gas with Z = 10~ 2 . For temperatures T < 10 4 K, we assume that 
ionization heating becomes important in such a way that d r T ~ 0. 

We integrated Eq. 18 above, as well as Eq. 9, in order to obtain both the velocity and 
temperature radial distributions. In Figure 4 we present the profiles of T and u/c s for two 
models with To = 10 5 K, L = lOOpc, NFW gravitational potential with po — 10~ 21 g cm -3 
and Rs = 200pc. The difference between the two is the initial energy density of waves that 
was set as (SV 2 ) Q /c 2 = 10 (solid lines) and 20 (dashed lines). It is obvious that the cooling 
and adiabatic expansion initially result in a temperature decrease while, on the other hand, 
the damping of the wave flux heats the gas. At lower radii the effect of the wave damping 
is negligible and the temperature gradient is negative for both cases. We found that the 
radiative cooling will dominate over the adiabatic expansion as main coolant of the plasma 
if the gas density is larger than n ~ 3cm -3 . 

As r — > L, there is an increase of the energy released by the waves for the plasma 
heating. In both models shown in Fig. 4 the wave heating becomes larger than the cooling 
processes and an increase in temperature is observed. For {5V 2 ) /c 2 = 20 the temperature at 
r ~ lOOpc becomes even larger than the initial setup temperature of To = 10 5 K. Depending 
on the initial wave energy density chosen it is possible to obtain temperatures as high as 
10 6 K. In both cases though, the temperature at the second critical point is T ~ 10 4 K. 
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Figure 4. Wind sonic Mach number u/c s ,o and gas temperature T as a function of the radius, for (SV 2 ) = 10c 2 s (solid) and 
20c 2 s (dashed), both with To = 10 5 K and L = lOOpc. 

Again we found no significant differences in terminal wind velocities, being Moo/c s ,o — 2.2 
and 2.4 for {SV 2 ) j<? s = 10 and 20, respectively. The reason is that most of the wave energy 
is deposited below the second critical point. However, notice that the Mach number is much 
different, M. oo ~ 7.5, because of the lower temperature at larger radii. 

Similar to what is observed in the isothermal cases, the fact that most of the wave energy 
density is deposited below the second critical position is relevant for the mass loss rate. We 
obtained M ~ 1.2 and 1.7 x 10~ 4 M Q yr _1 for {5V 2 } j<? s = 10 and 20, respectively. The mass 
loss rates obtained for the non-isothermal cases are more than 1 order of magnitude smaller 
than the obtained for isothermal winds with similar physical parameters, with T = 10 5 K. 
The obtained mass loss rates though are similar to the observed mass loss rates of colder 
isothermal winds (T ~ 10 4 K). This because the main acceleration of the non-isothermal 
winds occurred near to the second critical position, i.e. r > 400pc, where the temperature is 
already low, close to 10 4 K. We run several models of non-isothermal winds and all presented 
very low mass loss rates as well. If any relevant external heating source is included though, 
in the sense that the temperature does not diminishes substantially at r < r crit 2 , the mass 
loss rates would happen to be much larger, comparable to the isothermal models with larger 
temperatures. It is interesting that, despite of its increased complexity, the wind profiles 
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obtained from non-isothermal models are not relevantly different from those obtained with 
the simplified isothermal set of equations. 

4 DISCUSSION 

The results obtained from the semi-analytical modeling of galactic winds presented in this 
work show that dwarf galaxies are very likely to present steady winds with significant mass 
loss rates, but with relatively slow terminal speeds. Even in the case with low wave en- 
ergy density, i.e. basically pure thermal winds, the mass loss rates reach values as large as 
lO^nfooM^yr -1 , depending on the plasma temperature and the dark matter potential pro- 
file. This value can even be an order of magnitude larger if kinetic wave pressure is added 
at the bottom of the galactic wind. Therefore for the case of an initial total gas mass of 
M gas ~ 1O 7 M , as considered in the present work, the galaxy would loose the totality of its 
gas mass in timescales of order of 10 8 — 10 1 0yrs. This can be directly compared to the obser- 
vational estimates. Schwartz & Martin (2004), for example, detected outflowing gas from six 
dwarf galaxies with starbursts at very low velocities, typically Uqo ^ 40km s _1 ~ c s , as pre- 
dicted from our models. They also infered the outflow gas masses as M outflow ~ 10 6 — 1O 1O M , 
also in agreement with the total gas mass that such a system could loose during its evolution 
based on our model. Large mass losses from dSph have also been estimated for the Local 
Group objects based on chemical evolution models (Kirby et al. 2011). Slow galactic winds 
for dwarf galaxies have also been been detected by other studies, e.g. van Eymeren et al. 
(2010). 

Obviously, the mass loss rate is increased once wave pressure is added at small radii. 
Kinetic wave pressure is typically related to SNe explosions. Few explosions of the massive 
stars from young stellar populations would result in a turbulent interstellar gas, or even 
in the direct blow of the ejecta out of the galaxy. Supernovae (SNe) are possibly the most 
important source of turbulence in the ISM of dwarf galaxies, which have been generically 
considered in this work as a preset energy density of acoustic waves. Let us consider Type 
II SNe as the most important energy sources in an early dwarf galaxy. In a stationary case 
the wave energy flux equals the total power released from the stellar explosions into kinetic 
motions: 



-E'wav — £-^SN 
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(19) 



being 7£ S n the supernova rate, e SN ~ 10 51 erg the total energy released by the SN, and £ the 
fraction of the total energy converted into kinetic motions of the interstellar medium. 

Notice that for £ = 0.01 and n = 10 3 cm" 3 it is possible to obtain (5V 2 ) ~ 10c;? with 
TZsn ~ 10~ 7 SN yr _1 . This is a very low SNe explosion rate, compared to the estimates for 
typical early evolutionary stages of dSph (TZsn ~ 10 _6 SN yr _1 ). Therefore, it is reasonable to 
assume the existence of a continuous and isotropic acoustic wave pressure, from supersonic 
turbulence, at the central regions of the galaxy. 

Many models have been proposed in the literature to account for the mass loss rates of 
SNe driven winds using simplified prescriptions (e.g. Holzer & Axford 1970, Larson 1974, 
Bradamante et al. 1998, Lanfranchi & Matteucci 2007, Cote et al. 2012, Sharma & Nath 
2012, and others), but none have presented a systematic study of quasi- stationary solutions 
as has been presented in this work. 

Compared to previous works, an important difference from our results is the dependence 
of the wind mass loss rate with the profile of dark matter. This result shows that simply 
comparing the volumetric integrations of SNe energy/mass feedback and gravitational bind- 
ing energy may not the best relation to determine if the wind will be launched or not. This 
because the radial dependence of the acceleration of gravity is strongly dependent on the 
dark matter profile. Actually, this parameter is the one to be compared to the gradient of 
the energy density released by SNe. Another basic result, which gives different predictions 
compared to previous works, is the possibility of stationary winds even in the presence of 
very low star formation rate. Large mass loss rates are also obtained from pure thermally 
driven winds, given that any internal or external source keeps the interstellar gas heated 
at minimum temperature T > 10 4 — 10 5 K. Therefore, the quenching of star formation may 
occur at very early stages of the galaxy evolution. Estimates predict a reionization of early 
Universe at z ~ 5 — 11 (e.g. Robertson et al. 2010) due to the formation of the first stars, 
which corresponds to a cosmic age of 400 — HOOMyrs in a standard ACDM Universe. If the 
heating of the galactic ISM occurs at these times, and considering the timescales for mass 
loss from our models, we may estimate a cut-off in the star formation at t ~ 1.5Gyrs. 

Other key timescales for the galaxy evolution, such as the star formation timescale and 
the gas infall timescale, also play a major role on the time evolution of the equilibrium gas 
fraction since both timescales evolve strongly with redshift (see Khochfar & Silk 2011). The 
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accretion rate related to the cold gas infall is often related to the dark matter potential 
as M acc ~ /&Mdm, being /& ~ 0.165 the baryon mass fraction (Dekel et al. 2009). For the 
dwarf galaxies observed at z — > it is expected that their dark matter haloes did not change 
their masses much since their formation, except for merging processes. Larger galaxies, on 
the other hand, would evolve for longer times with an approximate growth rate given by 
M DM (z) ~ 0.033/ 6 ~ 1 M 1 1 o 15 (l + 2) 2 - 25 M yr" 1 (Neistein et al. 2006), where M 10 represents 
the dark matter halo mass in units of 10 1 solar masses. At z — 9, for a 10 9 solar masses 
dark matter halo, one then obtain r acc (z = 9) ~ /^Mdm^cc ~ 3 x 10 7 yr. The accretion is 
then supposed to occur much faster than the gas loss timescales given above. The typical 
timescales for star formation, on the other hand, is estimated as r SF ~ f b MDy[SFR~ l , 
being SFR the star formation rate. At z — 9, for a 10 9 solar masses dark matter halo, the 
SFR ~ 10~ 3 - lO^Meyr" 1 (see Weisz et al. 2011), which leads to r SF ^ 10 10 yrs. The two 
later timescales show that the intergalactic baryons must fall into the dwarf dark matter 
haloes, but the timescales needed for complete conversion into stars are much larger than 
the typical mass loss timescale. As a result, only few percent of the gas mass forms stars, 
and most of it is ejected out of the galaxy. Cases of galaxies with star formation bursts at 
2—t-O could be explained by later infall of gas, subsequent to a long period of inactivity. 

4.1 Dark matter and the time dependence of mass loss rates in dSphs 

Numerical simulations of the dynamical evolution of the matter in the early Universe predict 
the formation of a large number of low mass dark matter structures, which are supposed to 
host the dwarf galaxies. These simulations also predict a cusped dark matter profile in the 
center of these objects (e.g. Navarro, Frenk & White 1996). However, recent observational 
studies have pointed to less steep (cored) distributions instead (Burkert 1995; van den Bosch 
et al. 2000; de Blok & Bosma 2002; Kleyna et al. 2003; Simon et al. 2005; Walker et al. 2009; 
Governato et al. 2010; Oh et al. 2011; Del Popolo 2012; Jardel & Gebhardt 2012). 

The discrepancy between the cusped and core profiles is possibly related to the evolution 
of these galaxies. Since these low mass structures are typically satellites of massive galaxies, 
models predict that the interaction between their gravitational potentials stirs up the initially 
central-cusped dark matter distribution, resulting in cored profiles (Mayer et al. 2001a, 
2001b; Klimentowski et al. 2007, 2009; Kazantzidis et al. 2011; Lokas, Karantzidis & Mayer 
2012). Also, at later evolutionary stages of these galaxies the baryonic matter may interfere 
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in the dynamics of the dark matter. Possibly, stellar dynamics also result in core-like dark 
matter density profiles (Governato et al. 2010; de Souza et al. 2011; Pontzen & Governato 
2012; Governato et al. 2012). 

The timescales for changes in the dark matter distribution to occur range from r ch ~ 
hundreds of Myrs to Gyrs. These can be directly compared to the timescale for the galactic 
wind of a NFW profile to remove most of the interstellar gas r win d ~ M gas /MNF\v- If r ch > 
T win d, the sharp dark matter potential would drive a faster mass loss process that removes 
most the gas from the galaxy, and star formation is quenched until another gas infall from 
the intergalactic medium occurs (infall regulated star formation events). On the other hand 
if T ch < r wind , the gravitational potential evolves fast enough for the galactic mass loss to 
slow down. The interstellar gas would be only partially removed, and star formation would 
occur continuously. 

5 CONCLUSIONS 

Galactic winds play an important role on the mass feedback of the intergalactic medium, 
particularly on its chemical evolution, as well as in the galactic star formation history. Dwarf 
galaxies are known to present relatively more massive dark matter haloes, and systematically 
lower metallicities, compared to the more massive galaxies. The presence of strong galactic 
winds during the early stages of their evolution could be responsible for quenching star 
formation and help reduce the baryon fraction of their total masses. Theoretical modeling 
of dwarf galaxy winds typically considers supernovae (SNe) as main, or even only, sources 
of energy and momentum. 

In this work we apply a Parker-like wind model to the physical properties of dwarf galax- 
ies in order to study the role of the thermal pressure on their wind acceleration. We studied 
different dark matter potentials, temperature distributions and SNe driven turbulence. De- 
spite of the general picture of chaotic and stochastic galactic winds (e.g. SNe driven winds) 
it is possible to drive quasi- stationary thermally driven flows. Also, we showed that the sta- 
tionary solution presents two critical points, which will fully determine the wind mass loss 
rates and wind terminal velocities. 

As main results, we found that: 

• physical quasi-stationary solutions of winds in dwarf galaxies are possible if the inter- 
stellar gas is hot (T > 10 4 K), depending on the gravitational potential; 
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• dwarf galaxies are likely to present galactic winds at very early evolutionary stages, 
even when presenting very low star formation, since thermally driven winds showed to be 
quite efficient; 

• cusped dark matter profiles, e.g. the Navarro- Frenk- White profile, result in larger mass 
loss rates, compared to cored or less steep profiles. 

The models also showed that increased turbulent motions at the bottom of the winds, 
possibly associated to star formation activity, are responsible for an increase in mass loss 
rate. This result is in agreement to the simplified analytical models that compare the binding 
gravitational energy to the SNe energy thermalized at the ISM, and with numerical simula- 
tions of SNe driven galactic winds. However, the wind terminal velocities obtained showed 
that this parameter is insensitive to the SF activity. 

The models described in this work may serve as benchmarks for future numerical simu- 
lations, as well as basic wind mass loss rates and terminal velocity prescriptions for chemical 
and dynamical models of dwarf galaxies that use galactic outflows as a model input. As an 
example, closed box chemical evolution models include infall and mass loss rates in order to 
study the metal enrichment of the galaxies taking into account losses of metals to the IGM 
due to galactic winds. Typically, this has been done by simple energy budget arguments. 
These models would greatly benefit from the semi-analytical predictions of this work. Also, 
full hydro dynamical numerical simulations could also be set in order to check the validity of 
a quasi-steady solution for the wave/thermally driven galactic wind. As basic parameters, 
the set of a dark matter potential, initial equilibrium distributions for gas temperature and 
density, and the sources for kinetic waves and heating, would be enough in order to obtain 
the time evolution of the gas flow, and to compare with our results. 
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